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Abstract 

Calculations are presented for the thermodynainic functions and phase 
separation boundaries of solid metallic hydrogen-helium alloys at temperatures 
between 0‘^K and 19,000'*K and at pressures between 15 and 90 megabars. Expressions 
for the band structure energy of a randomly disordered alloy {including third order 
in the electron-ion interaction) are derived and evaluated. Short- and long-range 
order are Included by the quasi-chemical method, and lattice dynamics in the virtual 
crystal harmonic approximation. We conclude that at temperatures below 4,000®K there 
is complete phase separation of hydrogen-helium alloys, and that a miscibility gap 
remains at the highest temperatures and pressures considered. Tbe relevance of these 
results to models of the deep interior of Jupiter is briefly discussed. 




1. INTnOOUCTION 

Knowledge of the phase diagram of hydrogen-helium alloys nt high 

pressures (4-40 megobars) Is of Importance In the study of the Interior of 
fi 1-3 

the giant plants. Phase separation of hydrogen and helium during the 

2 

cooling process may partly account for Jupiter's excess emission of energy, 
this piper presents a calculation of the thermodynamic functions and phase- 
separation boundaries of solid hydrogen-helium alloys at pressures between 
15 and 90 megnbars, and at temperatures between O'^K and 19,000^K. These 
metallic systems are also of Intrinsic interest, since the particles carry 
point charges, and the bare electron-electron, electron-ion, and lon-lon 
Interactions are given exactljr by Coulomb’s law. 

3 

Xhe calculations reported here supplement earlier results of Stevenson 

on hydrogen-helium phase separation in the liquid phase. Pro sent estimates 

4 

of the melting curves of these materials and of the temperature in the deep 

3 

Interior of Jupiter Indicate that both hydrogen and helium may well be 

liquid in the planet’s interior, at temperatures far below 19,000®K. However, 

since the uncertainties In the calculated melting temperatures are quite 
5 

large , a solid-solid phase separation calculation remains of particular interest. 
The phenomenon of solid-solid phase separation in alloys is, of course, 

6 

not limited to the hydrogen-helium system, but is known to occur in many alloys. 

For example, Li and Mg (both simple metals) form solid alloys at all concentrations 
except in the range of about 707o-85% Mg, where there exists a miscibility 
gap. An alloy formed in this concentration range will separate into two phases 
of different concentrations. It is noteworthy that the miscibility gap is still 
present at temperatures Just below melting. The hydrogen-helirm alloy xs, how- 
ever, different from many other alloys (such as Li and Mg) in one Important respect. 


2 
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Whereas the difforonco between the Mg and Li oloctron-lon Interactions 

(pseudopotontials) is small, hydrogen and helium have olectron»ion interactions 

of very different strengths, and this difference is expected to play an 

inporthnt role in the thermodynamic pn^erties of their alloys. 

In Sec. II we discuss the general approach taken in formulating the 

Helmholtz free energy P for hydrogen, helium, and their alloys. The static 

Internal energy E is calculated in Sec. Ill for any given configuration of 
0 

hydrogen and helium (confined, however, to an underlying lattice), and is 
subsequently evaluated for a randomly - disordered configuration. Contributions 
to F arising from long - and short-range order are treated in Bee. IV, and 
the free energy associated with lattice dynamics in Sec. V. In Sec. VI we 
present the equations of state and the Gibbs free energy G per ion of hydrogen- 
helium alloys. Writing G as a function of its natural variables (pressure p, 
temperature T , and the relative concentration by number of helium c), we 
compute iCf which is defined by: 

AD = G(p,T,c) - [c G(p,T,l) + (1-c) G(p,T,0)] . (1) 

Fran AD we determine the curves describing solid-solid phase separation. 

II. HELMHOI.TZ FREE ENERGY 

For a system of volume 0, the free energy F can be written as 

F(T,n,c) := Fg(T,n,c) + Fy(T,Q,c) , (2) 

where F^(T,n,c) is the static free energy, and F^(T,0,c) the vibrational 
free energy. In principle, F can be calculated from the partition function Z, 
which is the sum of e ^ over all degrees of freedom, electronic and ionic, 
and in particular over all configurations of hydrogen and helium on the assumed 
underlying lattice. (Here g = and E is the total energy.) It is useful 

to introduce the following notation: Let (a) denote the ensemble average 

s 

of the variable A for a static lattice, The electronic degrees of freedom and 
the configurational degrees of freedom remain summed over in obtaining (a) • 


! 
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We use the symbol ^ to Indicate the ensemble average of A for a static 

lattice in which the configurations summed over are restricted to bo randomly 
disordered. Wo can now write F^(T, fj,c) of Eq, (2) as 

F^(T, n,c) = <1> - T <S> , (3) 

8 5 S 

where S is the entropy. We may also write Eq. <2) as 

•" =<'>.,» 4”. - k- "v] • 

where is the vibrational free energy of a randomly disordered alloy. 

W* will ignore the last term in Eq . (4), and in Sec. V calculate only 

F^ . Hie validity of this approximation will be discussed in the final section. 

Hie neglect of the term [F^ - F°] , and the separatl<»i of the static free energy 

as shown in Eq. (4), are motivated by the fact that those temperatures for 

which hydrogen-helium alloys actually do form are sufficiently high as to favor 

such rand can disorder. (This point will be argued more fully in Sec. VIO Thus 

we expect that at these temperatures (e) will be the major contribution to 

& I o 

<e) . Note that the second term of Eq. (4) includes the configurational entropy, 

s 

as well as corrections to the static energy due to correlations of the positions 
of hydrogen and helium on the lattice. 

III. STATIC ENERGY 

In this section we calculate <e) by writing a general expression for 

s , o 

E , the static energy of any configuration of the Ions, and then computing its 
s — - 

average over randomly disordered configurations. Hie approach is to consider 

an alloy as consisting of hydrogen and helium icms, located on a lattice, and 

immersed in a responding electron gas of compensating density. The ion-ion 

electron-electron, and electron-ion interactions are all given by Coulomb's 

law. Hie (divergent) long-wavelength limits of these interactions sum to 

7 

zero, and are eliminated from the starting Hamiltonian. One can then write 

E as 

s 

E = E^*^l E„ + E, . 

s Mb 


( 5 ) 


5 


^ !toro is the energy (per Ion) of n homogeneous Intornctlng electron gns (in the 

presence of a positive, uniform background charge), the Madelung energy E,. 

is the electrostatic energy of the point icms (in the presence of a uniform 

negative b.jckground charge), and £|^ is the energy due to the electrons' 

response to the non- uniform component of the total ionic potential Vf By 

treating V as relatively weak, Ej^, which is known as the band structure energy, 

can be calculated by perturbation theory. What we are describing is conven- 

8 

tional pseudopotential theory, applied to a system for which the electron-ion 

interaction is known exactly. 11118 approach has been used extensively in the 

7 9 

context of metallic hydrogen, * and is an important element in the alloy 
calculation of Ref. 3. 

In Eq. (5), is given by 

* 

where Z is the average ionic charge in units of e (e > 0) . Since Z „ » 2 

HE 

and Z_ = 1, Z - cZ,_ + (Z-c)Z = 1 + c. Note that r is the usual dimensionless 
H HE H S 

electron spacing parameter; 


(r a_)^ = 


(7) 


where a is the first Bohr radius. Since N is the number of ions (in D), NZ 
o 

is the corresponding number of electrons. Ihe first two terms (6) are 

the kinetic and exchange energies. Ihe last term is the correlation energy, and 

is only known approximately. We have used the approximation due to Nozieres and 
10 

Pines, which is expected to be quite satisfactory in the r range considered 

s 

here (r .s.. 1). Note that is independent of both the configuration of 

s 

hydrogen and helium ions on the underlying lattice, and of the lattice Itself. 

Since we are interested in temperatures much less than the Fermi temperature 

T„ = 5.82 X 10® K"" , (8) 

F 2 

r 

11 ® 

the electron system is taken to be in its ground state. 
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O 

Tho soco.id term in Eq . (5) is the Modolung energy, and is given by ; 


K„ = 1 
** 2T)N 


r r 

(iJ) !L 


4n 


ik.(R^-B^) 


ZiZj , 


( 9 ) 


where Z^ is the charge of tho Ion a- site i whoso position is given by R^. 

The prime cn tho sum over i and J denotes tho omission of the terms i » J. The 
prime on tho k-sum denotes tho omission of k » 0 . 


The Madelung energy is generally large and negative, and for a given f imily 

of structures often assumes its lowest value for the most symmetric structure, 

7 9 

Using perturbation theory ’ can be developed as a series in ascending 


orders of the electron~ion interaction: 


E - E<2> 4. E 4. 


with 


- i S I 'i B(i,) - '] , 

k, ?TT 


( 10 ) 


( 11 ) 


and 


„(3) in 7 . y. V(k.) V(k«) V(-k -l^)r 6 .1 .1 a ^ 

^ = 3 - J ^ eTq,) €(q^) cPva/ » 




where the primes denote the omission and ^ = - 1 ^ • 

In Eq.s (11) and (12), V(k) is given by: 

V 1 f j3 -ik.r,,, . 

V(k) a 7 = r d r e V(r) , 

~ w s ~ 


V(r) = - S Z, 




(13) 

(14) 


where V(r) is the total ionic potential as seen by the electrons, The 

restrictions on the sums in Eq,s (11) and (12) follow from the form of the 
7 

Hamiltonian. The vectors q are defined by q = k/2k_, where the Fermi wave- 
vector kp Is given by the relation 

kp = 3 tt^ Z*N . 

In Eq.s (11) and (12), €(q) is the zero-frequency limit of the dielectric function 

(3) 

of the homogeneous interacting electron gas, and H' ‘ is given in Eq. (C3) of 

0 


Ref, 7. Wo use Hartree atomic units in tho equations above, (and throughout 


( 12 ) 


the rest of the onoer) 


I 


7 

It Jb important to note that Eq. (11) is an exact roBult for 

b ’ 

**i[jG(q ) exact linear response of the number density 

of the homogeneous Interacting electron gas to an external potential, (in this 

case the potential due to the ions). In contrast, Eq. (12) is only approximate, 

as the corresponding second-order response function is not known exactly. Hie 

approximation used in £q.(12) corresponds to treating the electrons as independent 

particles moving in a self-consistent potential constructed from a Hartree potential 

and the external potential, provided e is taken to be the Lindhard dielectric 

function^ In the present calculations we have used the Ooldart-Vosko^^ 

modified form of th'^ Hubbard dielectric function, which includes effects due 

to exchange and correlation, and yields the correct q -* 0 limit. It is certainly 

preferable to use this form (rather than the Lindhard function) in e 5^\ but it 

b 

/O \ 

technically inconsistent to use it in E^ as written in Eq. (12). However, 

/ 3 \ 

these two dielectric functions yield values of E^ within 1% of each other, so 

b ’ 

that the effect on phase boundaries, which depend cm differences of free energies, 
is inconsequential. 

Although the hydrogen-helium alloys have been taken as metallic, the 

convergence of the perturbation series of Eq. (10) is not dependent on the 

existence of a metallic state, as discussed in Ref. 3. The point is that the 

perturbation series should be adequate as long as the one-electron band gaps 

are less than the band widths, which is the case for helium above 10 megabars. 

Since actual metallic conduction may only occur in helium at 70 megabars, this 

distinction is of considerable importance. (Hydrogen, on the other hand, is 
14 

expected to be metallic at pressures of a few megabars,) 

7 ^ 4 ^ 

Considerable progress has been made in calculating E^ , which however. 


I 
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wo do not incltido horo. For motallic hydrogen ' 1b smallor than ’ 

by rough' y a factor of ton , and It Includes the effects of the change i.r 

the chemical potential of the electron gas duo to the presence of the ions. 

To correctly calculate one must use finite-temperature perturbation 

theory, as discussed in Ref. 7. 

(2) (3) 

The terms E.., E' , and e' are valid as written for any configuration 
M b b — 

of hydrogen and helium, and cemtain contributions that depend both on the 
configuration and on the structure of the underlying lattice. More specifically, 
since the total potential V(jr) in Eq. (14) takes the form of a sum over sites, 

E|^ will contain the following classes of terms: 

(1) Structure Independent terms, that is, terms independent of configuration 

,( 2 ) 


and lattice structure. These arise from the terms 
(3) 

E^ ' In which all sites coincide. 


in e; 


and 


(11) Two-Body ■ or ion-ion terras, lliese ccmiprise the remaining terms in 

£|^ * , and the terms in ' for which only two site labels coincide. 


(ill) Ihree-Body, or ion-lon-ion terras. These arise from the terms in E 

' ^ W 

in which no site labels coincide. 

There are, of course, four-body terms and terms involving more than 

four ions, but these originate in higher orders of perturbation theory. 

Recognizing that E is also a sura over ion-ion terms, we con group together 
M 

contributions to E in Eq . (5) by the classes (l)-(lii) above, and obtain: 

a, r - 1 I 


(3) 


V 




3N 


'ijk "5k’ 5jj"5k^ + . . . 


(15) 


5i’5ij’5k 


Here the pi’inies denote restrictions forbidding the terms i = j in the two body term, 
and the terms i = k and j = k (but not i = j) in the throe body term. Note that 
the two- and throe- body potentials depend on density and on the identity of the ions 
ot sites i and j, (ns well as on the separations R. - R . ) . 


9 


All tormn in i: which nro indopondcnt oX conXlgurntion and lattice atructuro 

are included in 'Oio point about rewriting Eq, (5) aa in Eq, (15) ia 

aimply that by aumming over the electron degrooa of fi*eedora (at T^O“K), 
wo have been able to write aa a sum over (density-dependent) effective 
pair and three-body potentials, plus term dependent only on density. Hiis 
recasting of Eq. (5) is clearly valid for any configuration of hydrogen and 
helium ions, and is a conceptually useful alternative to Eq. (5). 

We now calculate the first term in Eq . (4), the static energy of a 
randomly disordered system; 


<E), 


= E 


(o) 


+ <E..>„ + + <E. + . . . 


(3)- 


(16) 


s,o ’^M'o ' b 'o '“b o 

To do this wo must first give the definition of randomly disordered . To this 


end we introduce the quantity p^: 

p^ = 1 , if site 1 is occupied by a helium ion, (17) 

p^ = 0, if Bite 1 is occupied by a hydrogen ion. 

From its definition, one can see that p^ obeys the following relations: 

(P^)” “ Pi» " “ • ' ' (18a) 

<P^> = e, (18b) 

where the average in Eq . ^18b) is over all configurations , Introducing the 


auxiliary variables d^: 


wo have 


Pi - 


<d^> = 0 . 


(19) 


( 20 ) 


Since p^ measures the probability that site i is occupied by a helium ion, 

d^ measures the deviation of that probability from its average value , In 

Eq. (9) for E,,, we write Z as 
M 1 


^i = Pi"FlE ^ 


( 21 ) 


10 


ThuB E will dourly Involve avorogoB of tho typo (p. p.) . In terms of 

M * J 

those correlation functions wo dof ino a randomly dlsorderod system as one 
for which tho n-th order correlation function factors according to : 

<P. P. • • • Pi )o “ W K ^Pf^o • • *^Pi = <Pi ) <Pi ) • • • <Pi ) r 
where 1^ 1^ ^ , . ^1^ 

•HiUB for the twO“Sito correlation function wo obtain* 

W**j^o* ^^i^ = c* if i j, <P^Pj>Q ® “ ^**1^0 " ® 

Since i s J ir excluded from £q. (9), wo immediately have: 

I I 

<E^> « Z*® y y 4n ^ik. (R. - R j ) 

^ ° ^ * 
i.J k 


1^)0 Mudolung energy of a randomly disordered alloy is that of a pure metal 

1.6 

of ionic charge Z , (corresponding to tho so-called "virtual crystal"), 
and can be calculated by well-known techniques. 

To calculate ^E. ) we must first use Eq.s (13) and (14) to write V(k) 

w O 


in terms of the variable p^: 


V(k) = 2.6 ”^~S.i (^p + (1-p ) 

~ i M Vk2 Cl ^ ^ ^ 'Trrr'' ' 

where R^ls again the position vector of site i. Introducing d^ via Eq. (19), 


we obtain: 


V(k) 




-ik.R 


~ (u (y + fli(y) , 


where 


ii (k) = - [< 


, ,, 4n Z„ 1 -4n Z 

+ » J = lETir ' 

kJi 0 k2 n 




and 


fU(k) = - ^ (z^jj, - Zjj) = 


Q 


( 22 ) 


(23) 


(24) 


(25) 


(26) 


(27a) 


(27b) 


u 


From Kq.s (11) nnd (17), wo flntlj 




Si 


( 28 ) 


and 


(v<)ti> v(-ii)>o =. I .'‘Si'Si L +‘Si-!l) <(s(iCj) + dja)%))(u(-!^) + djaK-kj)))o .(29) 


1 i 

From Eq. (20) wo ae« that the cross terms in Eq. (29) vanish. Usin^ the relation 

(30) 


I 


“ik.n. „ . 

e /w “ N Oj^ 

i 


where K is any vector oi the reciprocal lattice, we have: 


<V(k. ) V(-k, )>o ^ N 6,, „ U(k, 

~1 k, ,K ~1 


) U(-k, ) + ZZe^^’^i 


i J 


jCiU(k^)(!JJ(-k^)j< d^dj>o 
<31) 


in the Appendix, we prove the relationship: 

e-^iil’Sd <d.d > = N(c-c^) , 

id ^ 

Siibstitutln^ Eq.s (31) and (32) into (28), and using 


Lim i y -i i 3 r d^k 

n » or (2vr J ' 


(32) 


we have tho final second-order result: 

^ y I? E^) “O fe) " 


(K) 


whore Q = K/2k_ . In Eq . (33), the first term is just the second order band 
^ F 

17 * 

structure energy of a pure metal of ionic charge Z . Tliis virtual crystal result 
is not correct for a randomly-disordered system, because in Eq. (29) the terms 
in which tho sites i and j coincide must be handled separately. However, it is 
worth noting that tho virtual crystal result correctly gives the structural- dependence 
of since the second term in Eq, (33) is clearly Independent of botli the lattice 


structui’e and the configuration of hydrogen and helium on the lattice. 


12 


We hovo written in a form thot is quite similar to other 

fi 16 

exprosslonR in the literature, ' and have used a rather indirect method to 

do so. This method, however, avoids much of the confusicm that would otherwise 

,(3). 


arise in the calculation of which we now turn. 


Equation (12) for con be written in the following form; 


12 


■‘rill <v<ii> V(5a) V(!^)>„ XaC^.as.aj) \ 


♦ I I 

n r T V 


k k k 
^ 




(34) 


where the function is defined by direct comparison of Eq.s (34)ond (12), 

However, we shall never need the explicit exprostsion for , but only its 

symmetry properties. Hie form of the function H^^^(-q, ,q_) in Eq. (12) 

guarantees that is symmetric with respect to the interchange of any two 
7.12 

arguments. ' Using Eq.s (27) and (30), wo have 

.3 


<v(|^)v<i^)v<k,)> » N I k 6v ,k V ,k 5%) B(kc) Bcfe 


) 


+ 11 6 , 


^ u Sg (lCg)£U (1^) 

^ N C (k^)S2%,lCi> S) (kg) AO (|^> 

♦ S3<ki.kB.ka) flj(ki> ai(!^) aJCkg). 

.*hero wo have defined 

e ^!!a*54 <d^^d >j, , 


i J 


and 


r-«' 


S_(k, ,k„,kl = 'Lh 'L e • 5m e ‘ 2n <d,d ) 

^ ^ i m n ^ J ' 

These functions are shown in the Appendix to be 


(35) 


(36a) 


(36b) 


®2^~l'5a^ “ ^ ^ +l^,K ^ ’ 


(37a) 


I 


)3 


and 

‘>3<Si.Ss.!!a> - » \ +k,+k„K • 

Substituting Eq.s (35), (36) and (37) intti Eq. (34), and making use 
of the Bymmotry of wo obtain; 

f _ I 


(37b) 


^ [»^ ‘k^.J, ‘k3.k^ 


) 


^ ^ 








(C-3C® + 2c^ ^u(k^) m(i^) ^jj(kg)] )<a (2i-Sa'33> 
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(38) 


The first term in Eq. (38) is the third order bend structure energy of the 
virtual crystal. As befo.e, there are corrections to the virtual crystal result 
which have their origin in the coincidence of sites in Eq. (34), However, 
now the corrections are structure-dependent . To see this more explicitly, we 


recast in terms of the function of Eq. (12). By using the symmetry 

properties of with respect to Interchange of arguments (see Ref. 7), we 

can rewrite Eq. (38) as: 


.(3) 


<xf >>0 = - 8f 


9Tr\9n/ 


<^HE- V'l 




Si -22 ^®^Sl“S2^^”® ^^S-rS2^ 


e(q) * Q e(<3> 


r ) ii‘^> 

q-Q e(q-Q)/ s (q, 


Q) 


+ ^(c-3c^^2c^) 


^^HE"V Jd\ Jd^q2Uje(q^) * ‘ lSl-a2 ®<a.r32^^”s^^ % -520 , 


( 39 ) 


As before, Q - K/2k^, and the prime in the double sum means we omit 

Q s= 0, Q_ = 0, and Q = Q . Since the second term in Eq. (39) involves a 

sum over the reciprocal lattice, it is clearly structure-dependent, Eq. (39) 

(3 ) 

is our final result for (e^ >o . 


14 


The polynomials in c that opponr In Eq.s (33) and (39), (tho basic results 

of this section), are cumulant polynomials P^(c), familiar from the theory of 

19 

electron states in the tight-binding model of randomly disordered alloys. 

They are defined by the generating function 

«• 

it P (c) X = jM (1 - c + ce ) , (40) 

s=l ■ 

which gives, 

P, (c) - c, 

2 

P_(c) = c - c , 

^ 2 3 

Pg(c) = c - 3c + 2o , . . . . . (41) 

The cumulonts arise in both problems for the same reason, namely that the 
decoupling of the correlation functions, illustrated in Eq. (22), does not 
hold when two or more sites coincide, This point has been stressed previously in 
Ref.s 20 and 21. 

IV. LONG AND SHORT RANGE ORDER 

We now turn to the second term in Eq. (4), namely F - (e) . In the 

s s, o 

previous section we have summed over the electronic degrees of freedom to 

obtain an effective Hamiltonian for the ions (Eq. (15)). The static partition 

— E/k T 

functiem (and nence the static free energy) can be obtained by summing e B 
over all (static) configurations of hydrogen and helium ions on the underlying 
lattice. To carry out this sum, we need a convenient language with which to 
describe the configurations. At high temperatures, this is acheived through 
the use of the correlation functions * etc. introduced 

in See. III. In general, a helium ion may be more likely to have a hydrogen 
ion as a nearest neighbor than another helium ion (or vice-versa) , but the 
probability (at high temperatures) of a very distant neighbor of the helium 
ion being another helium ion will depend only on the mean concentration of 


helium. 


15 


Tho corrolntlon functions (p^Pj) otc. arc Idonlly sultod to describe such 

short-range order, for we expect tho quantity ^iP^Pj) “ (Pj^)(Pj* to become 

very small us and become increasingly well separated. On the other hand, 

at very low temperatures, and particularly for stoichiometric compositions, 

the alloy, if it forms at all, is expected to take U' an almost ccmpletoly ordered 

state. (For example, if c- 0.5, the alloy may have the CsCl structure at T - 0**K.) 

It is clearly inappropriate to attempt to describe this situation with the 

correlation functions of tho type ^Pj^Pj)r since (Pj^Pj) “ ^Pj^)(Pj) is expected * 

to be infinitely long-ranged . Instead, it is convenient to introduce the notion 

24 25 

of long-range order, * which for the exam).le quoted above would be defined 
by the number of helium ions on '* *’ieht sites", i.e. the number of He ions 
on the "helium ion" sublattice. Hie point is, of course, that this number is 
1.00 at T = 0'’K. It also approaches rather abruptly the disordered value 
of 0.5 at the critical temperature above which there is no long-range 

order . 

Thus, any theory used to calculate P„ - (e) must be capable of describing 

* s I o 

these two very different types of behavior at low and high temperatures. More 
specifically, at low temperatures we have; 

him (P_ - <E> ^) = , (42) 

Where is the energy difference between the completely ordered phase and 

24 

its completely random counterpart. At extremely high temperatures we have 

Lim (F„ - <E>^ ) « -T <S> = k_T \ c ike + (1-c) i^l-c)l , (43) 

.j. _ s s,o s,o B L J 

where the expression on the RHS of Eq. (43) is s. ly the negative of the 
entropy of a randomly disordered alloy, weighted by the temperature. 

The first step in formulating such a theory is drastically to sirapllfly 


JSq. (15), and replace it by a nearest-neighbor model, viz: 


’/m ^HE-HE ^}IE-H ^ 


h-h] ' 
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whoro the sum Is ovor nearest neighbors only, and the pair interactions 

*JI-H ^ chosen to satisfy Eq.s (42) and (43). Note 

that since we are computing only the difference between energies, tne structure- 

independent term in Eq. (15) may be neglected. The appeal of the simple form 

in Eq. (44) is that It allows an exact mapping of the problem onto the anti- 
. 26,27 

ferromagnetic Islng model. In addition, the Hamiltonian of Eq. (44) 

24 

has received a great deal of attention as a model Hamiltonian of an alloy. 

Since wo only need keep terms dependent on configuration, it is easy to show 
that the pair interacticHis do not enter separately, but only in the standard 
combination, 


'' ~ V-H Vh **■ ' 

28 


(45) 


where v is assumed to be negative. 

The energy difference as calculated from the Hamiltonian of Eq. (44), 
is proportic«al to -v, with the proportionality constant depending on the 
(stoichiometric) composition and the assumed underlying lattice. It is there- 
fore compelling to choose v so that the energy difference ^ between ordered 

29 

and disordered alloy will be the true static energy difference, as calculated 
by the methods of Sec, III, l.e. with no restrictions to nearest neighbors. 
Providing our methods of solving tne model problem defined by Eq.{44) satisfies 
the limit in Eq.(43), the resulting function F(T, Q, c) - (e) will then 

5 p O 

cxliibit both the correct high and low temperature behavior. 

Such a method of solution of the model problem is provided by the quasi- 

2 5 30 

chemical approximation. ’ Ihe basic idea of the method is to treat clusters 
of ions as independent units, subject only to the conservation of the number 
of each t3’pe of ion consistent with a given long-range order, Tlie probability 


of a cluster having a certain configuration of hydrogen and helium ions is then 
simply given by the standard Boltzmann factor. 


X7 


If tho cluster Is chosen to bo the whole crystal, the result is exact. For 

smaller clusters, <in particular lor n few atoms), error Is introduced because 

the fact that a given site may bo port of two (or more) clusters is Ignored 

In assigning a probability that tho site Is occupied by (say) a helium atom. 

Nevertheless, tho method does toko into account correlation effects in o manner 

reminiscent of classical liquid theory. The free energy can be written down as 

a functlcm of temperature and long-range order only, and is to be minimized with 

respect to the latter. Tho quasl-chcmical approximation is thus able approximately 

to describe both long- and short-range order within one context. 
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The approximatiem is related to more accurate methods in that it is 

32 

the first of a hierarchy of approximation which can be substantially developed, 
although the calculations become extremely involved. It is most readily applied 
in tho following cases: 

(1) c s 0.5, where the underlying lattice is bcc, and the assumed ordered 
state is the CsCl structure. 

(il) c = 0.75 (or c= 0.25), where tho underlying lattice is fee, and 

the assumed ordered state is the Cu^Au structure. 

The method correctly predicts that for c = 0«25 alloys (ii), the order-disorder 

25 

transition is of first-order, that is, the long-range order drops discon- 

tinuously to tero at T . It also correctly predicts that the transition for 

c 

alloys of type (i) is of second order, with the long-range order vanishing 

continuously at T^ . The existence of short-range order above the transition 

temperature, and hence a configurational contribution to the specific heat, 

33 

is also described by the method, but the details of the experimental specific 

25 27 

heats are reproduced only qualitatively, ’ 


18 


When compared to more accurate solutions of the Islng model, the quasi- 

34 35 

chemical method's prediction of T is only very roughly correct. ' However, 

0 

calculation shows that in the very low temperature region the quantity F -(e) 

s s ,o 

for c s 0.5 agrees fairly well with the low temperature ising model series 

. 26 
. expansion. 

We have used the quasi-chemical approximation to calculate F - (e) 

s s,o 

for c = 0.25, c = 0.50, and c = 0.75 alloys by using the solutions corresponding 
to the categories (i) and (ii) above. The parameter v was chosen to yield 
the true static energy difference between ordered and disordered phases, 
as previously described. However, the assumed structures for the ordered 
and disordered phases in the calculation of ££ were chosen by criteria to be 
explained in Sec.s V and VI, and wore not consistent with the structures for 
which the quasi-chemical method was evaluated (see (i) and (ii) above). In 
addition, the c<xitributions of lottice vibrations and the third order band structure 


energy to ££ were neglected 
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These approximations are expected to have a serious 


effect near T , but should make little difference well above or below T . 

c’ c 
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Since 


A E is a function of r , we have constructed an approximate form for P (T,r ,c) -(e) 

B 888^0 

which has the correct high and low temperature limits. We have not assumed 
that the order-disorder transition occurs at constant volume, for the actual 
behavior of the alloys is determined in Sec. VI from the Gibbs energy G computed 
at constant pressure and temperature . 


V. LATTICE VIBRATIONS 

To calculate the contribution to the free energy of the lattice vibrations 

we first assume that the alloy is randomly disordered. The "phonon" spectrum 

of the random alloy is then calcv'ated by replacing each ion with one of charge 

Z .. and mass M The values of Z and M are chosen so that the long- 

efi — eff exx eff 

• 38 39 

wavcleugtli limit of the phonon spectrum is given correctly, ’ 
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4 


This is readily iteen to require 


“eff • “ “'“he 


( 46 ) 


■nd 


"eff - " -'"he • 

Itie force conatanta for an alloy of arbitrary configuration are defined 
(to second order In the electron-ion interaction) from Eq. (15): 

(1 ^ J) (47) 




’a ’8 «lj* 


R - R -R, 


There are three x^qiea of force constants (corresponding to hydrogen-hydrogen, 


hydrogen-helium, and helium-helium pairs), and from Eq.s (ll)-(H) these are 

♦ "e <5* - ‘as ‘!U - 5.J*. 


and 


~ H-HE 

‘ as <5i - !s> 

~ HE-HE 

aB 


'h V ‘as - 5s> 


(4S> 


‘ <5, - »,) ‘ as ‘Si ■ Sj’ 


Here 4 .(R) depends on r and may be written as; 
a p ~ s 


‘as‘S> ■ ’a ’S 


4n 1 -Ik.R 

k7 Too ® 


(49) 


In terms of force constants, Eq. (46) Is equivalent to the replacement of the 
three types of force constants with a particular type of "average" force 
constant . 

llie concept of phonons In disordered systems in general, and more 

specifically the use of average masses and force constants, has met with some 

success when applied to alloys whose constituent elements have similar masses 

20,40,41 

or force constants. Clearly the masses and force constants of pure 

hydrogen and helium are not close to each other, but some Justification for 

the replacement of an alloy by an "equivalent" pure system Is given by the 

20 

"virtual crystal approximation" for the phonon Green function. More specifically. 
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if we start with a pure syatom of point ions having mass and charge given by 

Eq. (46), and introduce the difference between the physical charges and masses 

and the ''average" ones as a perturbation,^^ then within this approximation 

the perturbation causes no change in the phonon Green function. 

We have evaluated the dynamical matrix of the pure system defined by 

Eq. (46) in the adiobatic and harmonic approximations, with the electron-lon 

interaction taken Into account up to socond order. This has been repeated for 

a variety of crystal structures and concentrations, including pure hydrogen 

43 

and helium. From the phonon frequencies, we calculate the vibrational free 
energy ^^2 


F° = kjjT 2^ ^;(^2sinhCi f h 


i)- 


(50) 




Whore 3 = Oi(<i;j) is the phonon frequency of wave number ^ and branch 

index J, and the sum I*; over the first Brillouin zone, This zone sum was 
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carried out using the special-point technique ' with a modest number (^ 10) 
of special points. 

Note that by using the harmonic approximation, the frequencies appearing 
in Eq. (50) depend on r but not mi temperature. In order for them to acquire 

S 

a temperature dependence, a more sophisticated approximation, such as the self- 

43 

consistent phonon theory, would be needed. However, some thermal expansion 

is included by using the harmonic frequencies, for the contribution of to 

the pressure is not negligible (see Fig.s (1) and (2)). 

The calculation of the phonon frequencies of the (randomly disordered) 

alloys and of hydrogen and helium was used as a guide in the choice of the 

lattice structure chosen for the calculatiMis of Sec. 111. The point is that 

these Coulomb systems (in the virtual crystal phonon approximation) are very 

29 

often harmonically unstable, as discussed by Deck and Straus. 


(By an instability, wo rofor to tho occurronco of imaginary phonon f roquonclos , ) 

Tbo lattice structures used in tho calculations of Sec. Ill, as described in 
detail in Sec. VI, wore chosen to give real frequencies. It should be noted, 
however, that tho relationship between instabilities in tho virtual crystal 
approximation and those in the real (randomly disordered) alloy is not clear. 

We shall assess the effect of our approximate treatment of tho phonons on tho 
phase boundaries in Sec. VI. 

VI. RESULTS AND DISCUSSKSI 
A. Choice of Lattice Structures. 

Here wo discuss tho lattice structures chosen to calculate the various 

contributions to Eq. (4). Ttis static energy differences between lattices are 

7 9 

in general very small, ’ especially when compared to the energy in the phonon 

system. (However, those energy differences may not be small compared to tho 

differences in phonon energies between lattices.) Hiis raises the question of 

whether these materials can ever solidify in the conventional sense, it should be noted 

that the energy differences are also not necessarily small when compared to the 

difference tG of the Gibbs energies between the alloy and the pure hydrogen and helium 

systems, as Fig, 3 Illustrates. An extensive search in Bravais lattice space 

for the structure of lowest energy (ac carried out in Ref. 9) is not feasible 

for this problem: we limited ourselves to the consideration of the bcc, fee, and 

hep (with variable c/a ratio) lattices in the calculations of ( E ) and F ° 

s,o V 

in Eq. (4). (Simple cubic lattices are harmonically quite unstable for these 
systems . ) 

For the randomly disordered alloys (and for pure hydrogen and helium), 

* 

either fee or bcc proved to be stable for all Z except in the range 
★ 

1.20 ^ Z 1,30, and the stable lattice was chosen for the calculations. 

* 

At Z = 1.25, hep (with c/a = 1.7) was stable, and this structure was therefore 
chosen in the concentration range near Z* = 1.25, *1110 lattices used to compute 
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. 1 ' 


(e) and F ° arc «ummarl?.od in Table I. The absence of an entry for 
s,o V 

a particular contribution to the energy indicates that the value of that con- 
tribution was obtained by interpolation from its values at other concentrations. 

Note that was calculated for fee, not hep, in the region 1.10 Ci Z ^1.35, 

It is not expected that this procedure will cause any significant error in the 

phase separation curves, in addition, the designated phases for Z - 1.00 and 

46 

1.2S are harmonically unstable at low densities (corresponding to prossures of 

less than 20 and 30 megabars, respectively). Previous calculaticnii^^ show that 

such Instabilities will only occur at much higher values of r (lower pressures) 

8 

when the phonon spectrum is calculated in the self-consistent harmonic theory. 

ITtus we adopted the procedure of extrapolating the phonon frequencies to lower 
density to calculate at low pressure. 

We now discuss the lattice structure of the ordered alloys used in 
calculating F - (E) by the methods described in Sec. IV, Hjc energy 
difference ^ E between ordered and randomly disordered states was calculnted 
for c = 0.25, 0.50, and 0.75. (For pure hydrogen and helium, A E, as well 
as F^ - <E) ^ clearly vanishes.) For the alloy of ® = 0.50, we have considered 
two types of lattices: 

(i) Simple tetragonal (st), with o basis of one helium and one hydrogen 

, 47 

ion, situated so that when c/n = 1.0, this lattice has the CsCl structure, 
(il) Face-centered tetragonal (fet), with a basis of one helium and one 

hydrogen ion, situated so that when c/a = 1.0, this lattice has the NnCl 
structure . 

As the fet lattice proved unstable for a wide range of c/a values, we jsed 
the st lattice at c/a *=1.0, where it is stable. 

We considered two structures for the ordered c = 0.25 (c = 0.75) alloys: 

(i) Simple tetragonal (st) lattice of helium (hydrogen) ions with a four- 
point basis, Tl:e helium (hydrogen) ion resides at the lattice point, 

and three hydrogen (helium) ions sit at the face centers, if all the 
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iuri8 wtiro idonticnl, tho Inttlco would bo fnco-contorod totmironal . 

(nils Is tiio Konornlistntion of tho Cu^Au structure to c/a ^ 1»00) 

(11) Oody-contorod totrngonul (bet) lattice of liollum (hydrogen) Ions with 

a four-point basis, llio hollum (hydrogen) Ion rosldes at tho li{> co 

point, and throe hydrogen (helium) Ions sit at the face-centers and udge 

midpoints. If all the Ions wore Identical, the lattice would be simnle 
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tetragonal, with half the original lattice constant. 

Of these two structures, the st lattice with c/a = 0.7 proved, for c » 0.75, 

to have the lowest static energy (to second order In the electron-ion Interoction) . 

Since this structure Is harmonically stable, the difference between its static 

energy and that of tho corresponding disordered alloy of Table 1 (bcc) was sot equal 

to AE, as required in the application of the quasi-chomlcal theory of Sec. 'V. For 

c = 0.25, neither of the two structures are harmonically stable (over a wide 
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range of c/a values). Ttiis may be a dynamic Indication of immisclbllity 
at T = O^K, or alternatively It may Indicate that these structures are energetically 
quite far from the structure an ordered alloy actually assumes. Of these 
two structures, tho bet lattice with c/a = 1.0 has the lowest static energy 
for r > 0.920 (p <28.9 megnbars at T = O^'K) , but tho st lattice with c/a « 1.0 
has tho lowest energy for r^ < 0,920. Tho static energy differences between 
those structures and the corresponding random alloy (hep) were used for A E 
in the calculation of Sec. IV . Tn Table II we present tho critical temperature 
as a function of pressure for tho order-disorder transition, as calculated from 
Sec. IV. 

In order to determine how serious an error was made in neglecting lattice 

vibrations in the computation of A E, we computed for the CsCl-structure 

alloy at T = 0"K and r^ = 0.99. nie result is within 7% (0.001 a.u. per ion) 

of the correspcMidlng random alloy (bcc) result. The difference is small, even 

On the scale of A G. This also shows that our neglect of the term [F - F °] 

V V 

in Eq. (4) is quite justified. 
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D. Phase Separation# 

The equations of state of pure hydrof^en and helium are presented in 
Pig. 8 1 and 2. For hydrogen, at T « 0'*K, they agree well with Caron's results 
(see Ref. 29). 

Under conditions of constant temperature and pressure, the free energy 
to bo minimized is the Gibbs free energy G; 

G(p,T,c) = F(p,T,c) + P Oq , (51) 

where p is the presvure and the volume per ion. Stability of mixed phases 
is determined by A G: 

AC » G(p,T,c) - [c C(p,T,l) 4- (1 - c) G <p,T,0)3 . (52) 

Here c » 1 refers to pure helium and c s 0 to pure hydrogen. In order for there 

to bo any mixing, AG must bo negutivo. A miscibility gap occurs when A G 

is negative but the system can lower its Gibbs energy by separating into a helium- 

rich phase and a hydrogen-rich phase,®® lliis is onstratod in Fig. 3, whore 

wo present typical results for A3(p,T,c> at fixed p and T. At any concentration 

between c = c^^ and c ss c^ the system can lower its Gibbs energy by separating 

into a helium-rich phase at c ss c. and a hydrogen-rich phase at c «= c„, with 

the relative amounts of the two phases being given by number conservation. 

For such a partially separated system, the Gibbs function is given by the dashed 

line in Fig. 3. TJjo error bars in Fig, 3 refer to the estimated computational 
51 

error, not the error duo to the various physical approximations made. We have also 
shown typical static energy differences (to second order) between lattice structures 
in Fig. 3, from which the sensitivity of the phase boundaries to lattice structure 


can be estimated. 
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The phase separation curves themselves are presented In Figs. 4-B, 

Note that the tempera turos for which mixing occurs ore generally well above 

the ordor'^d .Border transition tomporaturos listed in Table, ii. Thus, as 

we have mentioned, the details of this transition are not very important in the 

calculation of the phase boundaries, Tlio uncertainties in A G are the cause 

of the uncertainties in the phase boundaries, indicated by the cross-hatched regions. 

regions. The most striking features of the results ore (i) arc the 

persistence of a large miscibility gap at the highest temperatures and pressures, 

52 

and (il) the largo temperatures necessary for ony mixing to occur. The 

occurrence of large mixing temporntures is not dependent upon the approxim ^tions 
we have used to take into account short-range order and lattice vibrations, 
although the precise values of the mixing temperatures clearly are. The prediction 
of complete phase separation at temperatures below some temperature T reflects 

nj 

the largo positive values of A G for the static alloys (A G ~ kr,T ), In contrast, 

u m 

the large miscibility gap is primarily duo to the "plnring" of the phase boundary 
near c s 0,25, This is caused by the exceptionally low values of A G for 
c = 0.25 (see Fig, 3) at high temperatures, an effect for which the lattice 
dynamics is entirely responsible. 

The relatively low phonon frequencies predicted by the virtual crystal 
approximation for tho c = 0,25 randomly disordered alloys should be compared 
with the imaginary froq\>encios found for tho c = 0,21 l> ordered alloys. In both 
cases the alloy exhibits phonons whose f requoncics-squarod are low. This results, 
in one case, In a true instability, and in tho other case the low energy and 
high entropy resulting from these low frequencies greatly favor mixing. In respect 
of tho c s 0.25 alloys, it appears that tho treatment of tho lattice dynamics 
may bo quite crucial. 


I 
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A moro correct trontmont of the t!l«ordorotl nlloy (within the harmonic theory), 
and the application on the tempo rnturo-dopondcn t Bolf-consiatont (harmonic) 
phonon theory for oxnaplo, '^dy prociuco qualitative difforoncon in the phaao 
boundaries. OiiS such difference might be the disappearance of the miscibility 
gap at temperatures below 19,000’-'K. 

In conclusion, the calculation iredicts that until the temperature has 
reached a fairly high value, which will certainly depend upon pressure, there 
is essentially complete phase separation in solid alloys of metollic hydrogen 
and helium. This may bo regarded ns a fairly firm result, since it is not dependent 
in any crucial way upon the approximations used to compute 6 G. If hydrogen 
und helium are solid in some region of the interior of Jupiter, these conclusions 
have 0 direct bearing on any phase separation model of energy emission. 

We also predict a largo miscibility gop that persists to T = 19,000‘^K and 
p = 00 taegabnrs. However, this prediction depends upon the approximations we 
have used in treating the lattice dynamics of the alloys, and might well bo 
substantially modified by a wore detailed troatraent of the phonon spectrum. 

The third-order terms in the band-structure energy have little effect, tending 
to raise A G by only a email amount. Thus the approximate response function used 
in as well as the neglect of is not expected to have any 

important effect on the phase boundaries. The same is true of the use of the 
quasi-chemical approximation. 
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APPENDIX 


The cfilculatlon of ^E. (e Soc. II! requires the ovnluntlon 


of tho following n'.oragos! 


® I Z o“%‘24 o <d^dj> 

1 J 


d'o ' 


(A-1) 


and 




®3 ^ 


-ik, ,R. -ik_.R -ik.,R \ 

e ~jt o e '»n vd .d d / 

X tn n o 


(A-2) 


i m n 

We will freely make use of the definitio* s and properties of tho variables 


p^ and d^ as presented in Sec, III, Expressing d^ in terms of p^, wo hove 


^‘‘i'^J^o ” <^Pi “ - c)>^ « <PiPj>o " 


(A-3) 


Similarly, 


<“lVn>o = «'’i - - ‘=><P„ - o»^ 


2 

= <P«P„P„>„ -c <p„n„>„ - c (p,p„>„ - C <p,p„>^ + 3c ~ 


(A-4) 


Xmno mno xno imo 

Note that if i m ^ n in Eq. ('-4), Eq. (22) guarantees that the average 
will vanish. If only two of the sites are equal, wo use Eq. (18) and again 
the average vanishes . Thus 

m ® + 2c^) . (A-5) 

f m n o X,m m,n 

Substituting Eq.s (A-3) and (A-5) into (A-1) and (A-2), and using Eq. (30), 


and 


=2<!li’Sa> ■ "V,!.„K > • 

<vL 'MS ^ 


®3^~l'Sa'}S3^ ” ,K + 2c^) , 

—•1 • J. — o ~ 


(A-6) 


(A-7) 


Where K is any vector of the reciprocal lattice. 
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of Eq. (46) is not exact, even in the long wavelength limit. The resulting 
error in the compressibility is of order 10%, (E. Stoll, P. Meier, and 

T. Schneider, II, Nuovo Cimento 23B, 90 (1974).) This discrepancy is 
also present in the case of pure hydrogen and helim. 

39. II. Beck and D. Straus (see Rof. 29) define the "average mass" incorrectly. 
However, since the mass of a pure system enters the dynamical matrix only 
as a multiplicative prefactor, none of their results are affected. 

40. W. A. Karaitakahara and B. N. Brockhouse, Phys. Rev. 1200 (1974). Note 

that the "average" force constants used in this reference do not correspond 
to the average defined by Eq.s (46) and (48), 

41. E. C. SvensBon, B. N. Brockhouse, and J. M. Rose, Sol. State Comm. 

245 (1965); S. C. Ng and B. N. Brockhouse, Sol. State Comm. 79 (1967). 

42. This procedure is necessary to keep r constant. 

43. P. Choquard, TTie Anharmonic Crystal (Benjamin, N.Y., 1971). 

44. A. Baldereschi, Phys. Rec. B7, 5212 (1973); D. J. Chadi and M. L. Cohen, 

Phys. Rev. B8, 5747 (1973). 

45. D. M. Straus and N. W. Ashcroft, Phys. Rev. BK, 448 (1976). 

4b . The type of "Kohn anomaly" instability shown by these two substances is 
discussed in Beck and Straus (see Ref. 29). The self-consistent phonon 
theory might well stabilize these substances at low density. 
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47. In tho context of cubic lattices, c/a is the ratio of the distanco 

botwoon equivalent planes to the distance between equivalent ions in a plane. 

48. P. Dyson, Ann. Phys. 1 (1971). 

49. Instabilities occur at long wavelength for both structures. 

50. Wo describe tho criterion for global instability. 

51. The largo error bars at higher temperatures and low concentrations of 
helium are largely due to the (estimated) error in only using a few special 
(hep) points to calculate P° for c = 0.25. The fractional error AF°/F*^ 

is usually less than 5%, but can be large, on the scale of (F^ 

for c = 0,25 In Fig. 3 is of order 0.1 a.u. per ion.) 

52. These features should be contrasted with the phase separation curves of 


Ref. 3. 
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FIGURE CAPTIONS 


Figure Equation of state of metallic hydrogen. 

Figure 2: Equation of state of helium. 

Figure 3: Typical results for AG vs. c. Ilie dashed line determines 
the phaso-separated region, (c^ < c < c^^). The dotted line 
shows another possibility for the phase-separated region 
consistent with the error bars. Typical static energy 
differences between lattices of randomly-disordered alloys 
are also shown. (FCT refers to face-centered tetragonal.) 

Figure 4: Phase separation curve at 15 megabars, x is the relative 

concentration (by number) of helium. The cross-hatched regions 
show the uncertainty in the phase separation boundary. 


Figure 

5; 

Phase 

separation 

curve 

at 

21 

megabars , 

Figure 

6: 

Phase 

separation 

curve 

at 

30 

megabnrs . 

Figure 

7: 

Phase 

separation 

curve 

at 

60 

megabars . 

Figure 

8; 

Phase 

separation 

curve 

at 

90 

megabars . 


TABLE I 


Lattices used in Computations for Randomly Disordered Alloys, and for Pure Hydrogen and Helium 


z 

1.00 

1.05 

1.10 

1.15 

1.20 

1.25 

1.30 

1.35 

1.40 

1.45 

1.50 


fee 

fee 

hep* 

hep 

hep 

hep 

hep 

hep 

fee 

bcc 

bcc 


fee 

fee 

fee 

fee 

fee 

fee 

— 

— 


— 

bcc 


"cc 

— 

— 

— 

— 

hep 

— 

— 


— 

bcc 

z 

1.50 

1.55 

1.60 

1.65 

1.70 

1.75 

1.80 

1.85 

1.90 

1.95 

2.00 


bcc 

bcc 

bcc 

bcc 

bcc 

bcc 

bcc 

bcc 

bcc 

bcc 

bcc 


bcc 

— 

— 

— 

— 

bcc 

— 

— 

— 

— 

bcc 

"v 

bcc 

— 

— 

— 

— 

bcc 

— 

— 

— 

— 

bcc 


*hcp refers to the hexagonal close-packed lattice with c/s “ 1.70. 
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TAOli: II 

3 

Order - dlHordcr critical temperature T (In units of 10 ”K) as a 

c 

function of pressure p (In units of megabars). (Pressures are approximate only). 


c 

T 

c 

■ 0.250 

P 

c ■ 

T 

c 

0.500 

P 

c 

T 

c 

- 0.750 

P 

5.06 

2.0 

3.45 

3.0 

0.79 

2.5 • 

4.82 

4.5 

4.40 

7.0 

1.21 

7.0 

4.65 

7.5 

5.63 

13.5 

1.70 

14.5 

4.45 

13.0 

6.67 

21 .0 

2.16 

23.5 

4.40 

20.5 

7.92 

34.0 

2.73 

39.0 

4.37 

31.0 

9.19 

50.0 

3.07 

49.5 

5.35 

47.5 

10.05 

63.5 

3.46 

64.0 

5.94 

59.5 

11 .03 

80.5 

3.89 

82.5 

6.61 

74.5 

12.10 

102.5 

4.21 

98.5 

7.35 

94.5 

12.68 

116.0 

4.47 

111.5 

7.90 

111.0 

13.31 

131.5 

4.75 

127.0 

8.33 

125.0 










p (Megabars) 



oocoi> 



p (Megabars) 



G/N (0.01 atomic units) 


I 








p « 21 Megabars 



p = 30 Megabars 














